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Summary 

Drag  caused  by  separated  flow  behind  the  hub  of  a 
helicopter  has  an  adverse  effect  on  aerodynamic  perfor¬ 
mance  of  the  aircraft.  To  determine  the  effect  of  sepa¬ 
rated  flow  on  a  configuration  used  extensively  for 
helicopter  aerodynamic  investigations,  an  experiment 
was  conducted  using  a  laser  velocimeter  to  measure 
velocities  in  the  wake  of  a  model  helicopter  hub  operat¬ 
ing  at  Mach-scaled  conditions  in  forward  flight. 

Velocity  measurements  were  taken  using  a  laser 
velocimeter  with  components  in  the  vertical  and  down¬ 
stream  directions.  Measurements  were  taken  at  13  sta¬ 
tions  downstream  from  the  rotor  hub.  At  each  station, 
measurements  were  taken  in  both  a  horizontal  and  verti¬ 
cal  row  of  locations.  These  measurements  were  analyzed 
for  harmonic  content  based  on  the  rotor  period  of  revolu¬ 
tion.  After  accounting  for  these  periodic  velocities,  the 
remaining  unsteady  velocities  were  treated  as  turbulence. 
Turbulence  intensity  distributions  are  presented.  Aver¬ 
age  turbulent  intensities  ranged  from  approximately 
2  percent  of  free  stream  to  over  15  percent  of  free  stream 
at  specific  locations  and  azimuths.  The  maximum  aver¬ 
age  value  of  turbulence  was  located  near  the  rear-facing 
region  of  the  fuselage. 

Introduction 

As  helicopters  demand  higher  forward  speeds,  ah 
sources  of  configuration  drag  must  be  minimized.  The 
adverse  effect  of  aerodynamic  forces,  principally  drag, 
produced  in  the  hub  region  of  helicopters  is  a  significant 
concern.  Data  obtained  during  helicopter  development 
show  that  the  percentage  of  drag  attributed  only  to  the 
hub  is  greater  than  22  percent  of  the  entire  aircraft  drag 
for  a  representative  configuration  (ref.  1).  Because  the 
hub  contributes  such  a  significant  portion  of  overall  vehi¬ 
cle  drag,  it  is  important  to  understand  and  accurately  pre¬ 
dict  the  sources  and  effects  of  hub  drag.  Another  adverse 
aerodynamic  effect  occurs  during  flight  conditions  when 
the  downstream  turbulent  flow  from  the  hub  interacts 
with  the  empennage  and  tail  rotor  in  standard  helicopter 
configurations.  This  effect  degrades  yaw  control  during 
these  flight  conditions. 

Modeling  the  helicopter  configuration  accurately  is 
important  for  preventing  significant  delays  during  vehi¬ 
cle  development.  The  complexity  used  in  modeling  the 
aerodynamics  of  helicopter  fuselages  varies  greatly. 
Parameters  may  be  determined  from  both  powered  and 
unpowered  scale  model  wind  tunnel  tests  and  from  anal¬ 
yses  ranging  greatly  in  complexity.  It  is  important  to 
understand  that  few  analyses  fully  account  for  the  differ¬ 
ences  in  fuselage  drag  caused  by  the  effect  of  the  rotor 
wake  on  the  fuselage,  an  effect  studied  during  wind  tun¬ 


nel  experiments  (ref.  2).  Currently,  only  experimental 
methods  using  powered  rotorcraft  fuselage  models  can 
give  insight  into  the  adverse  effects  of  the  hub  region 
described  in  the  previous  paragraph. 

Aerodynamic  investigations  were  conducted  to  mea¬ 
sure  the  velocity  fields  in  the  vicinity  of  a  model  helicop¬ 
ter  configuration  (ref.  3).  The  model  configuration 
consisted  of  a  generic  rotor  and  fuselage  representative 
of  a  wide  range  of  realistic  helicopters.  Comparisons  of 
the  measured  velocity  field  with  analytical  predictions 
were  made  (refs.  4  and  5)  to  help  understand  and  validate 
rotor  inflow  velocity  predictions.  The  effect  of  the  heli¬ 
copter  fuselage  on  these  inflow  velocity  measurements 
was  studied  (ref.  6).  However,  the  effects  of  the  region  of 
acknowledged  separation  behind  the  rotor  hub  were  not 
assessed. 

This  investigation  was  conducted  to  provide  and 
analyze  data  needed  to  form  the  basis  for  modeling  the 
region  behind  the  hub.  The  use  of  a  well-documented 
helicopter  model  configuration  used  for  previous  aerody¬ 
namic  configuration  studies  adds  to  the  value  of  this  spe¬ 
cific  experimental  data. 

Symbols 

Figure  1  shows  the  coordinate  system  used  in  this 
study. 

b  number  of  blades,  4 

T 

Cj  rotor  thrust  coefficient,  - - - - , 

0.0064  (nominal)  p%R~(QR)~ 

c  blade  chord,  2.61  in. 

F  Fourier  transform  of  once-per-revolution 

velocity 

F»  magnitude  of  harmonic  components  of 

u  velocity 

Fvr  magnitude  of  harmonic  components  of 

w  velocity 

Nu  number  of  points  in  each  azimuthal  bin  used 

for  computation  (represented  by  Nu  in  figures) 

Nw  measurement  count  for  vertical  velocity  mea¬ 

surement  (represented  by  Nw  in  figures) 

PMT  photomultiplier  tube 

R  rotor  radius,  2.82  ft 

r  radial  distance  from  hub  center,  ft 

T  thrust,  lbf 

Jl/2  (u'~  +  w'2) 

T  2D  turbulence  intensity,  - 

Moo 

u  velocity  in  x  direction,  ft/sec 

u  mean  velocity  in  x  direction,  ft/sec 


u  velocity  fluctuation  in  x  direction,  u  -  u, 

_  ft/sec 

2 

u'~  average  square  of  fluctuating  velocity  in 

x  direction 

free  stream  velocity,  ft/sec 
w  velocity  in  z  direction,  ft/sec 

w  mean  velocity  in  z  direction,  ft/sec 

w'  velocity  fluctuation  in  z  direction,  w  -  w, 

_  ft/sec 

2 

w'  average  square  of  fluctuating  velocity  in 

z  direction 

uy-  velocity  induced  in  forward  flight,  ft/sec 

w/,  velocity  induced  in  hover,  ft/sec 

x  downstream  distance  from  hub  center,  ft 

y  lateral  distance  from  hub  center,  ft 

z  vertical  distance  from  hub  center,  ft 

a  angle  of  attack  of  rotor  tip-path-plane  (posi¬ 

tive  nose  up),  deg 

p  density  of  air,  slug/ft 

a  standard  deviation 

\| /  rotor  azimuth  measured  from  downstream 

position  positive  counterclockwise  (as  viewed 
from  above),  deg 

Q.  rotational  speed  of  rotor  shaft,  221  rad/sec 

(nominal) 

Subscripts: 

/  forward  flight  quantity 

h  hover  quantity 

u  u  component  or  direction  quantity 

w  w  component  or  direction  quantity 

°o  free  stream  quantity 

Model  and  Apparatus 

Langley  14-  by  22-Foot  Subsonic  Tunnel 

The  Langley  14-  by  22-Foot  Subsonic  Tunnel 
(ref.  7)  is  a  closed-circuit  single-return  atmospheric  wind 
tunnel  (fig.  2).  The  unusual  test  requirements  associated 
with  vertical/short  takeoff  and  landing  (V/STOL)  and 
rotorcraft  aerodynamic  research  led  to  design  and  con¬ 
struction  of  this  tunnel  in  1970.  The  tunnel  has  a  test  sec¬ 
tion  that  can  be  operated  in  a  variety  of  configurations — 
closed,  slotted,  partially  open,  and  open.  The  closed  test 
section  is  14.5  ft  high  by  21.75  ft  wide  by  50  ft  long  with 
a  maximum  speed  of  about  338  ft/sec.  The  open  test  sec¬ 
tion  configuration,  which  has  a  maximum  speed  capabil¬ 
ity  of  about  270  ft/sec,  is  formed  by  raising  the  ceiling 


and  walls  to  create  a  floor-only  configuration.  For  this 
investigation  the  tunnel  was  operated  in  the  open  test  sec¬ 
tion  configuration,  where  a  nominal  turbulence  level  of 
0.25  percent  has  been  documented. 

The  tunnel  is  equipped  with  an  on-line  static  data 
reduction  system  that  can  display  computed  average 
aerodynamic  coefficients  with  interactions  and  wall 
interference  corrections  in  real  time.  A  significant  facil¬ 
ity  system  of  this  tunnel  is  the  two-component  laser 
velocimeter. 

Laser  Velocimeter 

The  laser  velocimeter  used  in  this  investigation  is 
designed  to  measure  the  components  of  velocity  in 
the  longitudinal  (downstream)  and  vertical  directions 
(ref.  8).  This  system  is  composed  of  four  subsystems: 
optics,  traversing  mechanism,  data  acquisition,  and 
seeding. 

The  optics  subsystem  (fig.  3)  operates  in  backscatter 
mode  and  at  a  relatively  high  power  to  overcome  the  long 
focal  lengths  needed  to  scan  the  width  of  the  tunnel  test 
section.  The  cylindrical  sample  volume  resulting  from 
the  beam  intersection  in  the  10-  to  20-ft  focal  length 
range  of  the  system  is  less  than  1  cm  in  length  with  a 
diameter  of  0.2  mm. 

The  traversing  mechanism  subsystem  includes  five 
axes  of  freedom — two  by  direct  translation  of  the  optics 
bench,  one  by  zoom  focus  of  the  sample  volume  across 
the  test  section,  and  two  by  tilting  of  the  final  steering 
mirror  as  the  laser  beams  exit  the  optics  package.  The 
traversing  mechanism  system  for  the  velocimeter  has 
nominal  downstream  travel  of  6  ft  and  vertical  travel  of 
12  ft.  This  travel  is  ideally  suited  for  scanning  the  entire 
flow  field  of  rotor  systems  less  than  6  ft  in  diameter. 

The  data  acquisition  subsystem  used  for  this  investi¬ 
gation  is  based  on  frequency  domain  processing  of  the 
signals  from  the  receive  optics  (ref.  9).  Velocity  infor¬ 
mation  from  the  frequency  domain  processor  and  azi¬ 
muthal  information  from  the  model  are  recorded  in  a 
buffer  interface  that  allows  for  later  azimuthal  correlation 
of  the  velocity  signals. 

The  principle  for  laser  velocimetry  is  based  on  anal¬ 
ysis  of  the  light  reflected  by  particles  that  move  with  the 
flow  through  a  specific  location.  The  seeding  system  in 
the  14-  by  22-Foot  Subsonic  Tunnel  (ref.  3)  inserts 
monodisperse  1 .7-pm-diametcr  polystyrene  particles  into 
the  portion  of  the  flow  where  the  laser  measurements  are 
being  made. 
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Rotor  Test  System 

The  rotor  test  system  used  for  the  experimental  data 
in  this  report  is  built  on  a  generic  test  system  (ref.  10) 
developed  at  the  14-  by  22-Foot  Subsonic  Tunnel.  This 
test  system,  the  2-Meter  Rotor  Test  System  (2MRTS), 
consists  of  a  29-hp  rotor  drive  system  with  collective  and 
cyclic  blade  pitch  controls,  a  four-blade  articulated  hub 
mounted  on  a  six-component  strain-gauge  balance,  and  a 
fuselage  skin  mounted  on  a  separate  similar  balance. 
These  two  six-component  strain-gauge  balances  provide 
independent  measurements  of  the  rotor  and  fuselage 
aerodynamic  loads.  Figure  4(a)  shows  the  test  configura¬ 
tion  installed  in  the  wind  tunnel  with  the  laser  velocime- 
ter  shown  in  the  background. 

The  shape  of  the  fuselage  is  designed  to  represent  a 
wide  range  of  helicopter  fuselages  without  being  specific 
to  any  particular  one.  The  fuselage  geometry  can  be 
described  by  a  set  of  equations  (refs.  6  and  1 1).  Develop¬ 
ment  of  computer  models  is  simplified  by  using  a  geom¬ 
etry  derived  from  equations.  Figure  4(b)  shows  the  shape 
of  the  fuselage  with  example  cross  sections. 


ation  when  gathering  laser  velocimeter  measurements. 
The  rotor  cyclic  controls  were  trimmed  to  eliminate 
once-per-revolution  flapping  of  the  rotor  blades.  Table  1 
lists  the  significant  flight  condition  parameters.  The 
peak-to-peak  deviation  and  standard  deviation  of  these 
parameters  were  determined  from  the  static  data  recorded 
at  each  measurement  and  are  listed  in  table  1  to  establish 
the  uniformity  of  the  flight  condition  from  one  velocity 
measurement  to  the  next.  Peak-to-peak  deviation  in  this 
context  refers  to  the  maximum  value  recorded  minus  the 
minimum  value  recorded. 

The  grid  of  points  where  measurements  were 
obtained  was  determined  from  consideration  of  mean 
rotor  downwash.  The  center  of  the  rotor  downwash  is 
expected  to  travel  on  a  trajectory  defined  by  the  forward 
speed  of  the  wind  tunnel  flow  and  the  average  momen¬ 
tum  imparted  to  the  flow  by  the  lifting  rotor.  From 
momentum  calculations  (ref.  1) 

4  3  114 

wf-2u00Wj ■  sin  a  +  m~  Wr-  wh  =  0  (1) 

where 


The  rotor  hub  used  in  this  study  is  a  four-bladed 
articulated  hub  (ref.  10).  The  hub  articulation  is  set  at 
2.0  in.  from  the  center  of  rotation  with  coincident  flap 
and  lag  hinges.  The  blade  grip  component  of  the  hub  has 
a  maximum  diameter  of  1.625  in.  and  extends  to  7.5  in. 
from  the  center  of  rotation.  The  airfoil  section  of  the 
blade  is  fully  developed  at  8.5  in.  from  the  center  of  rota¬ 
tion.  The  center  section  of  the  hub  has  a  maximum 
height  of  2.125  in.  from  the  plane  of  blade  articulation. 
The  pitch  rod  attaches  to  the  blade  grip  at  a  radius  of 
2.0  in.  from  the  center  of  rotation  and  is  attached  to  the 
swash  plate  at  a  radius  of  1 .4  in.  with  a  vertical  offset  of 
approximately  4  in.  Figure  4(c)  shows  the  exposed 
geometry  of  the  hub. 

The  rotor  system  tested  consists  of  four  rectangular 
blades  with  a  constant  NACA  0012  airfoil  section.  The 
blades  have  an  equivalent  linear  twist  of  -8.0°  from  the 
center  of  rotation  to  the  tip.  The  chord  of  2.61  in.  and 
radius  of  33.88  in.  give  the  system  a  solidity  ( bchiR )  of 
0.098. 

Test  and  Procedure 

The  velocity  measurements  were  taken  at  a  flight 
condition  similar  to  that  reported  for  a  previous  laser 
measurement  program  (ref.  12).  The  advance  ratio  was 
set  to  0.23.  A  nondimensional  thrust  coefficient  of 
0.0064  was  set  with  a  tip-path  angle  of  -3.0°  relative  to 
the  free  stream.  The  rotor  tip  speed  was  set  to 
624.0  ft/sec.  This  tip  speed  has  proven  to  be  benign  in 
terms  of  rotor  system  life  and  is  an  important  consider¬ 


Substituting  the  flight  conditions  of  this  investiga¬ 
tion,  the  vertical  component  of  forward  velocity,  Wf.  is 
8.61  ft/sec.  This  value  was  used  to  determine  a  uniform 
slope  for  the  grid  of  measurement  points.  Since  the  onset 
flow  is  approximately  144  ft/sec,  the  vertical  offset  of 
each  plane  of  measurements,  which  were  spaced  2  in. 
apart  in  the  downstream  direction,  is  approximately 
0.12  in.  This  vertical  offset  was  used  to  define  the  center 
of  each  downstream  grid.  There  were  nominally  15  mea¬ 
surement  locations  at  each  downstream  station,  1  at  the 
center  of  the  grid  defined  above,  4  spaced  above  the  cen¬ 
ter,  4  on  the  left,  4  on  the  right,  and  nominally  2  below. 
The  measurement  locations  are  given  in  table  2  and 
shown  in  figure  5. 

The  unsteady  flow  field  in  the  vicinity  of  a  helicopter 
rotor  in  forward  flight  is  expected  to  be  periodic  with 
blade  passage.  Measurements  taken  with  the  laser 
velocimeter  are  derived  from  the  random  distribution  of 
seed  particles  passing  through  the  sample  volume.  The 
resulting  random  interval  between  samples  requires  that 
a  special  technique  be  applied  to  the  data  to  determine 
the  periodic  content  of  the  velocity  measurement.  Since 
no  cross-stream  component  of  velocity  was  measured, 
the  swirl  induced  by  the  rotor  and  seen  in  the  lateral  com¬ 
ponent  of  velocity  was  not  measured. 

To  determine  the  periodic  unsteady  velocity  from 
laser  velocimeter  samples,  the  correlation  between  the 
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velocity  of  an  individual  sample  and  the  rotor  azimuth 
must  be  established.  Approximately  4000  samples  were 
taken  at  each  measurement  location.  The  laser  data 
buffer  stores  the  rotor  azimuth  with  each  valid  velocity 
measurement  from  the  frequency  domain  processor. 
Because  it  is  expected  that  the  rotor  will  produce  velocity 
that  is  periodic  with  blade  passage,  sorting  the  velocity 
measurements  by  the  azimuth  of  the  rotor  defines  the 
periodic  content.  For  this  analysis,  128  azimuthal  ranges 
or  “bins”  of  2.8°  were  used  for  the  sorting.  An  interac¬ 
tive  analysis  routine  for  validating  particle  velocity  dis¬ 
tributions  was  used  to  edit  the  measurement  data  prior  to 
entry  into  the  working  velocity  database.  This  routine 
allows  the  user  to  review  the  velocity  histogram  at  each 
measurement  point  and  eliminate  invalid  samples  such  as 
those  resulting  from  reflections  from  the  body.  After 
editing,  the  average  number  of  measurements  in  each 
azimuthal  u  velocity  bin  was  30.2;  the  average  number 
for  w  measurements  was  30.6. 

Repeat  measurements  were  taken  at  random  inter¬ 
vals  during  this  investigation.  Repeat  measurements 
were  also  taken  at  one  of  two  measurement  locations  at 
the  specified  test  condition  after  major  breaks  in  the  test 
process.  Data  resulting  from  these  repeat  measurements 
were  used  for  two  purposes:  (1)  confirmation  that  sys¬ 
tems  were  operating  properly  and  that  the  test  conditions 
were  maintained,  and  (2)  as  a  check  on  the  uncertainty  of 
the  measurement  technique.  When  discrepancies  were 
uncovered,  the  test  procedure  was  checked  to  determine 
the  source  of  the  discrepancies.  When  the  repeat  points 
did  not  discern  a  problem,  the  data  were  included  in  the 
Results  and  Discussion  section  to  demonstrate  the  gen¬ 
eral  repeatability  of  the  measurement  process. 

Measurement  Uncertainity 

Potential  sources  of  measurement  uncertainty  were 
considered.  Data  from  the  laser  velocimeter  suggest  sev¬ 
eral  effects  may  contribute  to  measurement  uncertainty 
and  are  given  as  follows: 


The  total  laser  velocimeter  error,  in  percent  of  velocity,  is 
estimated  as 


/  2  2 

Total  =  A/Bias“  +  Random”  =  1.33  percent 

This  estimate  presumes  the  actual  flow  accelerations  are 
no  greater  than  those  presented  in  reference  12. 

Analysis 

The  flow  in  the  region  downstream  of  a  helicopter 
hub  is  expected  to  have  two  principal  unsteady  compo¬ 
nents:  an  oscillating  component  with  period  equal  to 
rotor  blade  passage,  and  a  turbulent  or  aperiodic  compo¬ 
nent.  The  velocity  measurements  made  during  this 
experimental  program  were  analyzed  to  separate  these 
two  unsteady  components.  The  periodic  part  of  the 
velocity  was  determined  by  the  sorting  of  velocity  mea¬ 
surements  into  azimuthal  bins.  The  average  value  of 
velocities  in  these  bins  was  assumed  to  represent  the 
periodic  component  of  the  unsteady  velocity  at  the  mea¬ 
surement  location.  The  turbulence  at  a  measurement 
location  was  then  calculated  as  the  unsteady  deviation  of 
the  velocity  from  the  expected  periodic  value. 

Turbulent  intensity  was  computed  from  statistical 
information  gathered  during  the  laser  velocimeter  data 
reduction  process.  The  turbulent  intensity  was  computed 
as  a  function  of  azimuth  by  using  information  about  data 
in  each  azimuthal  bin.  If  the  turbulence  was  not  com¬ 
puted  as  a  function  of  azimuth,  the  expected  periodic 
variance  of  velocity  with  azimuth  would  influence  the 
level  of  computed  turbulent  intensity.  One  of  the  param¬ 
eters  of  the  data  reduction  process  is  the  standard  devia¬ 
tion  a  of  particles  in  each  azimuthal  bin  described  above. 
The  standard  deviation  a  is  computed  from 


Source  of 
uncertainty 

Bias, 

percent 

Random, 

percent 

Crossbeam  angle 

error 

±0.81 

Not  applicable 

Diverging  fringes 

Not  measured 

Not  measured 

Velocity  bias 

0 

0 

Bragg  bias 

0 

0 

Velocity  gradient 

Negligible 

Negligible 

Frequency  domain 
processor 

Not  applicable 

0.20 

Particle  lag 

±0.50 

Negligible 

Total  error 

±1.31 

0.2 

where  Nu  is  the  number  of  points  in  each  azimuthal  bin 
used  for  the  computation.  Since 


72 


u  = 


N„ 


(4) 


? 

then  the  term  u  is  computed  from  the  deviation  as 


,2  2 

u  =  <7 


V  Nu  , 


(5) 
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Turbulent  intensity  T  is  defined  (ref.  13)  by 


T  = 


i  /'rl  .  /2  ,2 

1/3  u  +  V’  +  w 


(6) 


For  two-dimensional  velocity  data  the  intensity  is  rede¬ 
fined  as 


T  = 


/2 


+  w 


u 


oo 


(7) 


From  this  follows  the  calculation  of  turbulent  intensity  T 
at  each  azimuthal  bin  from  the  statistical  data  computed 
from  the  velocity  measurements  in  each  bin  as  follows: 


T  = 


2  TV,  -  1 
T"  2 N 


+  cr 


_2  Nw  ~  1 
2  N... 


(8) 


The  value  of  mean  turbulence  at  every  measurement 
location  has  been  computed  as  the  average  of  the  turbu¬ 
lence  from  each  of  the  azimuthal  bins. 

Although  most  azimuthal  bins  had  sufficient  num¬ 
ber  of  particles  for  this  statistical  analysis,  a  few  mea¬ 
surement  locations  had  no  particles.  To  prevent  the 
assumption  that  no  particles  indicates  no  turbulence,  a 
linear  interpolation  of  the  turbulence  across  the  empty 
bins  was  used  to  prevent  skewing  the  mean  turbulent 
intensity  values. 


Results  and  Discussion 

Figure  6  shows  the  distribution  of  mean  turbulent 
intensity  in  the  lateral  and  vertical  directions  at  each 
downstream  station.  Figure  6  is  divided  to  show  15 
downstream  stations  (labled  a-o).  Each  downstream  sta¬ 
tion  in  figure  6  shows  two  parts:  the  left  part  of  the  figure 
is  the  lateral  distribution  of  mean  turbulent  intensity  and 
the  right  part  of  the  figure  shows  the  vertical  distribution 
of  mean  turbulent  intensity. 

Figure  6  shows  that  the  level  of  mean  turbulence  is 
clearly  greater  on  the  side  of  the  hub  where  the  rotor 
blades  are  retreating  relative  to  the  oncoming  flow.  On 
the  side  of  the  hub  where  the  blades  are  advancing,  the 
average  level  of  turbulence  is  less.  The  lateral  intensity 
levels  appear  to  reduce  slowly  as  the  flow  progresses 
downsteam  but  retain  the  relative  distribution  of  strength. 
Much  less  mean  turbulence  is  observed  above  the  rotor 
(z  >  0)  than  behind  or  below  the  rotor  (z  <  0).  The  value 
of  mean  turbulence  downstream  of  the  hub  center  (z  =  0) 


seems  to  decay  at  a  slower  rate  than  below  the  down¬ 
stream  center.  Beyond  24  in.  (r/R  =  0.73)  downstream, 
the  average  turbulence  directly  behind  the  hub  center  has 
a  larger  value  than  the  turbulence  values  below  it.  The 
significantly  lower  turbulence  values  seen  beyond 
y  -  5  in.  (r/R  =  0.15)  on  the  advancing  side  indicate  little 
or  no  separation  of  the  flow  at  the  root  of  the  blade  on  the 
advancing  side  at  this  flight  condition. 

Multiple  measurements  were  taken  at  two  specific 
locations.  (See  figs.  6(a)  and  6(i).)  In  the  location 
approximately  8  in.  above  the  hub  where  no  separation  is 
expected  (fig.  6(a)),  all  of  the  measurements  correlate 
well.  In  a  region  of  expected  separation  (fig.  6(i)) 
slightly  less  correlation  is  shown  although  the  values  are 
consistent.  This  demonstration  of  repeatability  in  esti¬ 
mated  turbulence  from  the  statistics  of  velocity  measure¬ 
ments  improves  confidence  in  the  method  of  obtaining 
turbulence  values. 

Figures  7-21  present  details  of  the  velocity  measure¬ 
ments  and  derived  turbulent  intensity  values  for  each 
measurement  location.  Each  figure  shows  seven  plots. 
The  top  two  plots  show  the  azimuthal  variation  of  the 
two  components  of  velocity.  The  following  two  plots 
(Fu  and  Fw)  are  the  spectral  amplitudes  of  the  velocity  as 
functions  of  the  order  number.  These  values  are  com¬ 
puted  by  performing  a  Fourier  transform  on  the  azimuth- 
dependent  velocity  and  computing  the  magnitude  at  each 
order  number.  The  order  ratio  represents  the  frequency 
ratio  of  the  periodic  velocity  signal  in  multiples  of  rotor 
fundamental  frequency;  for  example,  an  order  ratio  value 
of  2  indicates  a  2-per-rev  frequency  content  in  the  veloc¬ 
ity.  The  next  two  plots  show  the  number  of  independent 
measurements  taken  of  the  two  velocity  components  at 
each  azimuth.  The  final  plot  in  each  figure  is  the  com¬ 
puted  turbulence  intensity  value  as  a  function  of  rotor 
azimuth. 

From  the  velocity  and  azimuthal  information, 
several  details  of  the  flow  behind  the  center  of  the 
model  rotor  system  can  be  seen.  For  example,  in  fig¬ 
ures  7(a)-(e)  the  downstream  component  of  velocity 
demonstrates  an  increasing  level  of  twice  per  revolution. 
This  2-per-rev  character,  although  unexpected,  has  been 
observed  during  previous  tests.  This  characteristic  is 
observed  at  downstream  stations  (figs.  7-10,  (a)-(e)) 
until  the  measurements  at  16  in.  behind  the  hub.  At 
16  in.  (r/R  =  0.49)  behind  the  hub  the  blade  passage  (four 
per  revolution)  frequency  content  exceeds  the  2-per-rev 
content  at  locations  above  the  hub  centerline. 

A  possible  significance  of  4  per  rev  exceeding  2  per 
rev  is  that  the  magnitude  of  the  load  on  the  blade 
(four  per  revolution  in  the  fixed  system)  becomes  signifi¬ 
cant  at  16  in.  from  the  center  of  the  hub.  After  16  in. 
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downstream  of  the  hub  the  2-per-rev  character  shifts 
closer  to  the  fuselage  at  most  stations  (figs.  1 1-2 1(e),  (f), 
and  (g)).  There  are  also  some  measurement  locations 
where  the  existence  of  2-per-rev  occurs  without  an 
upstream/downstream  correlation,  as  seen  in  fig¬ 
ures  12(k),  13(i),  17(h)  and  (i),  and  18(i).  Although  the 
2-per-rev  seems  to  be  most  closely  related  to  flow  associ¬ 
ated  with  the  interaction  of  the  hub  flow  and  the  fuselage, 
no  definite  explanation  is  currently  available  to  com¬ 
pletely  explain  its  existence.  A  possible  explanation  is 
the  existence  of  geometry  on  the  hub  that  is  2  per  rev, 
meaning  two  opposed  flapping  potentiometers  and  a 
rotating  scissors  that  is  mirrored  about  the  shaft.  The  rel¬ 
ative  frontal  areas  of  these  geometric  features  are,  how¬ 
ever,  relatively  small. 

Beginning  with  figure  10,  the  flow  associated  with 
the  back-facing  region  on  the  fuselage  can  be  seen. 
Scales  for  the  downstream  component  of  velocity  at  the 
two  stations  closest  to  the  fuselage  have  been  adjusted  to 
show  the  reverse  flow  that  occurs  at  the  two  stations  clos¬ 
est  to  the  fuselage.  At  the  next  station  downstream,  the 
downstream  velocity  is  also  quite  low  (figs.  11(f)  and 
(g)).  Figures  6(d)— (i)  show  average  turbulence  values  as 
high  as  approximately  15  percent  and  azimuthal  peaks  of 
over  20  percent.  This  region  of  the  highest  measured 
levels  of  turbulence  extends  back  to  approximately  26  in. 
behind  the  hub  (77  percent  of  radius). 

Another  flow  feature  that  is  seen  only  in  the  detailed 
velocity  and  azimuthal  plots  is  the  periodic  content  com¬ 
parison  between  the  advancing  and  the  retreating  side  of 
the  rotor  hub.  For  example,  in  figure  13(h)  (9  in.  to  the 
right  of  the  hub  and  20  in.  downstream)  the  vertical  com¬ 
ponent  of  velocity  has  a  smooth  but  periodic  character. 
By  comparison,  figure  13(o)  (9  in.  to  the  left  of  the  hub 
but  also  20  in.  downstream)  also  shows  a  4-per- 
revolution  content  in  the  vertical  component  of  velocity, 
but  it  is  not  smooth.  The  turbulent  intensity  computed  on 
the  left  (retreating)  side  of  the  rotor  hub  experiences  a 
level  of  about  10  percent  of  free  stream  while  the  right 
(advancing)  side  experiences  only  about  a  2-percent  level 
of  free  stream. 

Figure  14(e)  includes  very  similar  measurements 
taken  at  two  different  times  during  the  flow  measure¬ 
ments.  Figure  15(g)  also  shows  four  separate  sets  of  very 
similar  measurements.  Although  the  number  of  samples 
collected  in  each  of  the  azimuth  bins  determines  the  qual¬ 
ity  of  the  statistical  data  derived  from  these  measure¬ 
ments,  the  observed  repeatability  of  measurements 
contributes  to  the  subjective  confidence  in  the  method. 

Close  interaction  between  the  rotor  blade  and  the 
measurement  location  is  seen  in  the  velocity  distribution 


at  stations  between  22  in.  and  36  in.  (r/R  =  0.67  to  1.09) 
downstream  of  the  hub  (figs.  14-21),  especially  at  the 
first  location  above  the  mean  downwash  point  (4  in. 
above,  figs.  14—2 1  (d)  in  all  cases).  At  azimuths  of  0°, 
90°,  180°,  and  270°  the  vertical  component  of  velocity 
shows  a  sudden  disruption.  The  character  of  the  vertical 
component  of  velocity  seen  in  figure  14(d)  indicates  a 
concentration  of  circulation  or  vorticity  in  the  flow  field 
close  to  the  measurement  location.  This  character  is 
believed  to  come  from  the  abrupt  change  in  the  bound 
circulation  at  the  rotor  blade  tip  as  it  passes  by  the  mea¬ 
surement  location.  This  is  further  confirmed  by  the  com¬ 
plete  lack  of  such  a  signature  in  the  downstream 
component  of  velocity.  This  probably  indicates  that  the 
circulation  is  aligned  with  the  flow  at  the  time  of  the 
measurement. 

In  figure  19(d)  for  station  32,  the  downstream  com¬ 
ponent  of  velocity  shows  a  perturbation  corresponding  to 
the  vertical  component  event.  The  blades  are  32.88  in. 
long  and  the  effect  of  the  end  of  the  blade  just  beyond  the 
measurement  location  of  32  in.  now  contributes  to  the 
downstream  component  of  velocity  as  the  blade  passes 
the  measurement  location.  This  character  confirms  the 
general  understanding  of  the  transfer  of  bound  circula¬ 
tion  to  trailed  circulation  manifest  in  the  tip  vortex.  This 
signature  is  also  seen  in  the  downstream  component  of 
velocity  in  figure  20(d)  for  station  34  and  figure  21(d)  for 
station  36. 

In  figure  14(d)  the  measurement  count  goes  to  zero 
at  three  azimuthal  locations.  In  figures  15-19  ((d)  for 
each  figure)  the  measurement  count  (Nw)  for  the  vertical 
velocity  measurement  goes  to  zero  at  four  azimuthal 
locations.  In  figures  18(d)  and  19(d)  both  the  vertical 
and  lateral  velocity  measurement  counts  go  to  zero  at 
four  azimuthal  locations.  This  lack  of  measurement  is 
believed  to  be  caused  by  a  blade  obscuring  or  penetrating 
the  measurement  volume  and  interrupting  the  laser 
velocimeter  signal. 

In  figure  12(h)  a  “spike”  is  seen  in  the  turbulent 
intensity  data  at  approximately  304°  of  azimuth.  Analy¬ 
sis  of  the  contributing  data  indicates  that  there  is  a  single 
“high”  data  element  in  the  deviation  of  velocity  in  the 
downstream  direction.  In  figure  12(h)  the  value  is 
22.7  fps  while  the  neighboring  data  are  on  the  order  of 
3  fps.  Since  the  number  of  particle  measurements  in  this 
azimuthal  bin  is  32,  the  reason  for  this  spike  is  unknown. 

Conclusions 

The  unsteady  flow  behind  the  hub  of  a  helicopter 
configuration  influences  the  overall  flow  characteristics 
of  the  helicopter.  The  unsteady  flow  behind  the  hub  of  a 
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research  configuration  was  measured.  This  configura¬ 
tion  has  been  documented  by  extensive  laser  velocimeter 
measurements  and  provides  an  excellent  data  set  for  ana¬ 
lytical  correlation.  The  measurement  of  the  unsteady 
character  of  the  flow  behind  the  hub  was  taken  to  assist 
in  determination  of  relative  effects  of  this  unsteady 
region. 

Flow  behind  the  hub  of  this  characteristic  helicopter 
configuration  was  measured  at  13  stations  downstream 
from  the  rotor  hub.  At  each  station  measurements  were 
made  in  both  a  horizontal  and  vertical  row  of  locations. 
The  downstream  stations  extend  beyond  the  radius  of  the 
rotor.  At  each  measurement  point  the  vertical  and  down¬ 
stream  components  of  velocity  were  measured  using  a 
laser  velocimeter.  Approximately  4000  samples  were 
taken  at  each  measurement  point.  For  each  sample,  the 
rotor  azimuth  was  recorded.  The  correlation  of  velocity 
measurements  with  rotor  azimuth  allowed  definition  of 
periodic  velocities  at  each  measurement  location. 

The  unsteady  flow  in  the  region  behind  the  hub  has 
been  analyzed  in  two  parts.  The  principal  part  of  the 
unsteady  flow  is  presumed  to  be  caused  by  the  periodic 
passage  of  the  rotor  blades.  The  secondary  part  of  the 
unsteady  flow  is  caused  by  turbulence  or  aperiodic 
flow.  An  unusual  characteristic  of  the  flow  is  the  twice- 
per-revolution  periodic  characteristic  found  in  specific 
regions  of  the  flow.  The  turbulent  flow  has  been  charac¬ 
terized  with  a  two-dimensional  turbulent  intensity  that  is 
independent  of  the  periodic  part  of  the  unsteady  flow. 

The  resulting  characterization  of  the  unsteady  flow 
behind  the  rotor  hub  shows  that  the  aperiodic  magnitude 
of  the  flow  is  greater  on  the  retreating  side  than  the 
advancing  side  of  the  rotor.  The  periodic  magnitudes  are 
shown  to  be  greater  on  the  advancing  side  of  the  rotor. 
There  is  a  region  of  the  largest  aperiodic  flow  behind  the 
rear-facing  slope  of  the  helicopter  fuselage  shape.  The 
magnitudes  of  the  aperiodic  flow  can  be  as  much  as 
15  percent  of  the  free  stream. 

The  velocity  field  associated  with  the  concentration 
of  circulation  on  the  rotor  blades  is  seen  in  the  azimuth- 
dependent  velocities  at  measurement  locations  close  to 
the  passage  of  the  rotor  blades.  Some  azimuthal  mea¬ 
surements  were  obscured  because  of  the  passage  of  the 
blades  at  these  locations.  Measurements  close  to  the  end 
of  the  rotor  blades  show  the  influence  of  the  abrupt 
change  in  the  bound  circulation  at  the  blade  tip  as  it 
passes  by  the  measurement  location. 


NASA  Langley  Research  Center 
Hampton,  VA  23681-0001 
September  10,  1996 
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Table  1 .  Flight  Condition  Parameters 


Parameter 

Average 

Peak  to  peak 
deviation 

Standard 

deviation 

Velocity,  ft/sec 

143.4 

1.8 

0.36 

Tip  path  attitude,  deg 

-2.7 

1.1 

0.28 

Density,  Slugs/ft3 

0.002330 

0.000043 

0.000008 

Hover  tip  speed,  ft/sec 

621.4 

3.2 

0.58 

Advance  ratio 

0.230 

0.003 

0.0006 

Collective  cyclic,  deg 

6.4 

0.2 

0.05 

Lateral  cyclic,  deg 

-0.8 

0.5 

0.16 

Longitudinal  cyclic,  deg 

2.8 

0.6 

0.22 

Lift,  lbf 

145.6 

5.2 

1.1 

Drag,  lbf 

0.73 

1.23 

0.273 

Thrust  coefficient 

0.0064 

0.0003 

0.00006 

Power,  hp 

7.6 

0.9 

0.24 

Rotation  rate,  rpm 

2110.8 

8.1 

1.42 

Unit  Reynolds  number,  per  ft 

889874 

29060 

5673 

Table  2.  Measurement  Locations 


x,  in. 

y,  in. 

z,  in. 

Figure 

8 

0 

8.518 

7(a) 

8 

0 

7.518 

7(b) 

8 

0 

5.518 

7(c) 

8 

0 

3.518 

7(d) 

8 

0 

-.482 

7(e) 

8 

0 

-4.482 

7(f) 

8 

9 

-.482 

7(g) 

8 

8 

-.482 

7(h) 

8 

6 

-.482 

7(i) 

8 

4 

-.482 

7(j) 

8 

-4 

-.482 

7(k) 

8 

-6 

-.482 

7(1) 

8 

-8 

-.482 

7(m) 

8 

-9 

-.482 

7(n) 

10 

0 

8.398 

8(a) 

10 

0 

7.398 

8(b) 

10 

0 

5.398 

8(c) 

10 

0 

3.398 

8(d) 

10 

0 

-.602 

8(e) 

10 

0 

-4.602 

8(f) 

10 

9 

-.602 

8(g) 

10 

8 

-.602 

8(h) 

10 

6 

-.602 

8(i) 

10 

4 

-.602 

8(j) 

10 

-4 

-.602 

8(k) 

10 

-6 

-.602 

8(1) 

10 

-8 

-.602 

8(m) 

10 

-9 

-.602 

8(n) 

12 

0 

8.278 

9(a) 

12 

0 

7.278 

9(b) 

12 

0 

5.278 

9(c) 

12 

0 

3.278 

9(d) 

12 

0 

-.722 

9(e) 

12 

0 

-4.722 

9(f) 

12 

9 

-.722 

9(g) 

12 

8 

-.722 

9(h) 

12 

6 

-.722 

9(i) 

12 

4 

-.722 

9(j) 

12 

-4 

-.722 

9(k) 

12 

-6 

-.722 

9(1) 

12 

-8 

-.722 

9(m) 

12 

-9 

-.722 

9(n) 

I 


Table  2.  Continued 


x,  in. 

y9  in. 

z,  in. 

20 

0 

7.796 

13(a) 

20 

0 

6.796 

20 

0 

4.796 

20 

0 

2.796 

20 

0 

-1.204 

13(e) 

20 

0 

-5.204 

13(f) 

20 

0 

-7.204 

20 

9 

-1.204 

20 

8 

-1.204 

20 

6 

-1.204 

20 

4 

-1.204 

20 

-4 

-1.204 

20 

-6 

-1.204 

13(m) 

20 

-8 

-1.204 

13(n) 

20 

-9 

-1.204 

13(o) 

22 

0 

7.676 

14(a) 

22 

0 

6.676 

14(b) 

22 

0 

4.676 

14(c) 

22 

0 

2.676 

14(d) 

22 

0 

-1.324 

14(e) 

22 

0 

-5.324 

14(f) 

22 

0 

-7.324 

14(g) 

22 

9 

-1.324 

22 

8 

-1.324 

22 

6 

-1.324 

22 

4 

-1.324 

22 

-4 

-1.324 

22 

-6 
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Figure  3.  Schematic  of  laser  velocimeter  optic  subsystem. 


(a)  Helicopter  configuration  in  wind  tunnel. 
Figure  4.  Model  configuration. 
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(c)  Hub  geometry.  Dimensions  in  inches. 
Figure  4.  Concluded. 
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Figure  5.  Velocity  measurement  locations. 
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Figure  7.  Concluded. 
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Figure  8.  Continued. 
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Figure  9.  Velocity  and  turbulence  at  station  12  in. 
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Figure  9.  Continued. 
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(g)  y  =  9  in.,  z  =  -0.722  in. 
Figure  9.  Continued. 
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Figure  9.  Continued. 
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Figure  10.  Continued. 
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Figure  11.  Velocity  and  turbulence  at  station 
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Figure  11.  Continued. 
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Figure  11.  Continued. 
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(i)  y  =  8  in.,  z  =  -0.963  in. 
Figure  11.  Continued. 
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90 


180 

Azimuth,  deg 


-50.0 


J _ 1 _ 1 _ 1 _ 1 _ I _ 1 _ L 


Fu 


Fw 


Nu 


Nw 


92 


10.0 


Orders 


Figure  11.  Continued. 
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Figure  12.  Continued. 
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Figure  12.  Continued. 
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(o)  y  =  - 9  in.,  z  =  -1.083  in. 
Figure  12.  Concluded. 
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Figure  13.  Velocity  and  turbulence  at  sta 
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Figure  13.  Continued. 
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Figure  13.  Continued. 
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Figure  13.  Continued. 
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Figure  13.  Concluded. 


Fu 


Fw 


Nu 


Nw 


T 


10.0 


1  23456789  10 

Orders 


Figure  14.  Velocity  and  turbulence  at 
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Figure  14.  Continued. 
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Figure  14.  Continued. 
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(g)  y  =  0  in.,  z  =  -7.324  in. 
Figure  14.  Continued. 
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Figure  14.  Continued. 
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Figure  14.  Continued. 
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Figure  15.  Velocity  and  turbulence  at  station 
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Figure  15.  Continued. 
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Figure  15.  Continued. 
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Figure  16.  Velocity  and  turbulence  at 
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Figure  17.  Velocity  and  turbulence  at  station 
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Figure  18.  Velocity  and  turbulence  at 
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Figure  19.  Velocity  and  turbulence  at 
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Figure  20.  Velocity  and  turbulence  at 
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Figure  20.  Continued. 
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Figure  20.  Continued. 


J _ I _ 1 _ I _ I _ I _ I _ L_ 


J _ I _ 1 _ 1 _ 1 _ 1 _ I _ L_ 


270 


360 


12  13  14  15  16 


223 


270 


Fu 


Fw 


Nu 


Nw 


Figure  20.  Continued. 


Fu 


Fw 


Nu 


Nw 


T 


10.0 


Orders 


Figure  20.  Continued. 


J _ I _ 1 _ I _ I _ 1 _ 1 _ L_ 


270 


360 


11  12  13  14  15  16 


J _ I _ 1 _ I _ I _ 1 _ I _ L 


360 


270 


Fu 


Azimuth,  deg 


Fw 


Nu 


Nw 


T 


10.0 


1  23456789  10  11 

Orders 


Figure  20.  Continued. 


J _ I _ I _ I _ I _ I _ 1 _ L. 


270 


360 


12  13  14  15  16 


J _ 1 _ I _ I _ 1 _ I _ ! _ L_ 


270 


360 


T 


50.0  p 
25.0  E- 
W,  fps  0  r- 


-25.0^ 
-50.0 


J _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ L 


90 


180 

Azimuth,  deg 


Fu 


10.0 


Fw 


Nu 


Nw 


T 


10.0 


1  23456789  10 

Orders 


Figure  21.  Velocity  and  turbulence  at 
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Figure  21.  Concluded. 
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